NCIO 1.0: A simple Fortran NetCDF interface

The NetCDF (Network Common Data Form) library has become an indispensable tool for data and model output management in geoscience. However for simple tasks, particularly in Fortran, the complexity of native NetCDF functionality can be cumbersome. The NCIO (NetCDF Input/Output) module has been designed as an interface to the NetCDF library with simplicity and ease of use in mind. While this implies that some NetCDF functionality is masked from the user, the subroutines provided here are adequate for basic serial reading and writing tasks of up to 6-D data arrays along with corresponding data attributes. The code is available online via a GitHub repository (http://www.github.com/alex-robinson/ncio), which includes an example program to illustrate the approach

Differential climate impacts for policy-relevant limits to global warming: The case of 1.5 °c and 2 °c

Robust appraisals of climate impacts at different levels of global-mean temperature increase are vital to guide assessments of dangerous anthropogenic interference with the climate system. The 2015 Paris Agreement includes a two-headed temperature goal: "holding the increase in the global average temperature to well below 2°C above pre-industrial levels and pursuing efforts to limit the temperature increase to 1.5°C". Despite the prominence of these two temperature limits, a comprehensive overview of the differences in climate impacts at these levels is still missing. Here we provide an assessment of key impacts of climate change at warming levels of 1.5°C and 2°C, including extreme weather events, water availability, agricultural yields, sea-level rise and risk of coral reef loss. Our results reveal substantial differences in impacts between a 1.5°C and 2°C warming that are highly relevant for the assessment of dangerous anthropogenic interference with the climate system. For heat-related extremes, the additional 0.5°C increase in global-mean temperature marks the difference between events at the upper limit of present-day natural variability and a new climate regime, particularly in tropical regions. Similarly, this warming difference is likely to be decisive for the future of tropical coral reefs. In a scenario with an end-of-century warming of 2°C, virtually all tropical coral reefs are projected to be at risk of severe degradation due to temperature-induced bleaching from 2050 onwards. This fraction is reduced to about 90% in 2050 and projected to decline to 70% by 2100 for a 1.5°C scenario. Analyses of precipitation-related impacts reveal distinct regional differences and hot-spots of change emerge. Regional reduction in median water availability for the Mediterranean is found to nearly double from 9% to 17% between 1.5°C and 2°C, and the projected lengthening of regional dry spells increases from 7 to 11%. Projections for agricultural yields differ between crop types as well as world regions. While some (in particular high-latitude) regions may benefit, tropical regions like West Africa, South-East Asia, as well as Central and northern South America are projected to face substantial local yield reductions, particularly for wheat and maize. Best estimate sea-level rise projections based on two illustrative scenarios indicate a 50cm rise by 2100 relative to year 2000-levels for a 2°C scenario, and about 10 cm lower levels for a 1.5°C scenario. In a 1.5°C scenario, the rate of sea-level rise in 2100 would be reduced by about 30% compared to a 2°C scenario. Our findings highlight the importance of regional differentiation to assess both future climate risks and different vulnerabilities to incremental increases in global-mean temperature. The article provides a consistent and comprehensive assessment of existing projections and a good basis for future work on refining our understanding of the difference between impacts at 1.5°C and 2°C warming

Long-Term climate change commitment and reversibility: An EMIC intercomparison

This is the final version of the article. Available from the American Meteorological Society via the DOI in this record.This paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to 1) quantify the climate change commitment of different radiative forcing trajectories and 2) explore the extent to which climate change is reversible on human time scales. All commitment simulations follow the four representative concentration pathways (RCPs) and their extensions to year 2300. MostEMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near-preindustrial values in most models for RCPs 2.6-6.0. The MOC weakening is more persistent for RCP8.5. Elimination of anthropogenic CO2 emissions after 2300 results in slowly decreasing atmospheric CO2 concentrations. At year 3000 atmospheric CO2 is still at more than half its year-2300 level in all EMICs forRCPs 4.5-8.5. Surface air temperature remains constant or decreases slightly and thermosteric sea level rise continues for centuries after elimination ofCO2 emissions in allEMICs.Restoration of atmosphericCO2 fromRCPto preindustrial levels over 100-1000 years requires large artificial removal of CO2 from the atmosphere and does not result in the simultaneous return to preindustrial climate conditions, as surface air temperature and sea level response exhibit a substantial time lag relative to atmospheric CO2. © 2013 American Meteorological Society.KZ and AJW acknowledge support from the National Science and Engineering Research Council (NSERC) Discovery Grant Program. AJW acknowledges support from NSERC's G8 Research Councils Initiative on Multilateral Research Funding Program. AVE and IIM were supported by the President of Russia Grant 5467.2012.5, by the Russian Foundation for Basic Research, and by the programs of the Russian Academy of Sciences. EC, TF, HG, and GPB acknowledge support from the Belgian Federal Science Policy Office. FJ, RS, and MS acknowledge support by the Swiss National Science Foundation and by the European Project CARBOCHANGE (Grant 264879), which received funding from the European Commission's Seventh Framework Programme (FP7/2007–2013). PBH and NRE acknowledge support from EU FP7 Grant ERMITAGE 265170

The Earth system model CLIMBER-X v1.0 – Part 2: The global carbon cycle

The carbon cycle component of the newly developed Earth system model of intermediate complexity CLIMBER-X is presented. The model represents the cycling of carbon through the atmosphere, vegetation, soils, seawater and marine sediments. Exchanges of carbon with geological reservoirs occur through sediment burial, rock weathering and volcanic degassing. The state-of-the-art HAMOCC6 model is employed to simulate ocean biogeochemistry and marine sediment processes. The land model PALADYN simulates the processes related to vegetation and soil carbon dynamics, including permafrost and peatlands. The dust cycle in the model allows for an interactive determination of the input of the micro-nutrient iron into the ocean. A rock weathering scheme is implemented in the model, with the weathering rate depending on lithology, runoff and soil temperature. CLIMBER-X includes a simple representation of the methane cycle, with explicitly modelled natural emissions from land and the assumption of a constant residence time of CH4 in the atmosphere. Carbon isotopes 13C and 14C are tracked through all model compartments and provide a useful diagnostic for model–data comparison. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing the historical evolution of atmospheric CO2 and CH4 but also the spatial distribution of carbon on land and the 3D structure of biogeochemical ocean tracers. The analysis of model performance is complemented by an assessment of carbon cycle feedbacks and model sensitivities compared to state-of-the-art Coupled Model Intercomparison Project Phase 6 (CMIP6) models. Enabling an interactive carbon cycle in CLIMBER-X results in a relatively minor slow-down of model computational performance by ∼ 20 % compared to a throughput of ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer in a climate-only model set-up. CLIMBER-X is therefore well suited to investigating the feedbacks between climate and the carbon cycle on temporal scales ranging from decades to &gt;100 000 years.</p

Climatic risks and impacts in South Asia: extremes of water scarcity and excess

This paper reviews the current knowledge of climatic risks and impacts in South Asia associated with anthropogenic warming levels of 1.5°C to 4°C above pre-industrial values in the 21st century. It is based on the World Bank Report “Turn Down the Heat, Climate Extremes, Regional Impacts and the Case for Resilience” (2013). Many of the climate change impacts in the region, which appear quite severe even with relatively modest warming of 1.5–2°C, pose significant hazards to development. For example, increased monsoon variability and loss or glacial meltwater will likely confront populations with ongoing and multiple challenges. The result is a significant risk to stable and reliable water resources for the region, with increases in peak flows potentially causing floods and dry season flow reductions threatening agriculture. Irrespective of the anticipated economic development and growth, climate projections indicate that large parts of South Asia’s growing population and especially the poor are likely to remain highly vulnerable to climate change

Long-term climate change commitment and reversibility: an EMIC intercomparison

This paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to: (i) quantify the climate change commitment of different radiative forcing trajectories, and (ii) explore the extent to which climate change is reversible on human timescales. All commitment simulations follow the four Representative Concentration Pathways (RCPs) and their extensions to 2300. Most EMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near pre-industrial values in most models for RCPs 2.6–6.0. The MOC weakening is more persistent for RCP 8.5. Elimination of anthropogenic CO2 emissions after 2300 results in slowly decreasing atmospheric CO2 concentrations. At year 3000 atmospheric CO2 is still at more than half its year-2300 level in all EMICs for RCPs 4.5–8.5. Surface air temperature remains constant or decreases slightly and thermosteric sea level rise continues for centuries after elimination of CO2 emissions in all EMICs. Restoration of atmospheric CO2 from RCP to pre-industrial levels over 100–1000 years requires large artificial removal of CO2 from the atmosphere and does not result in the simultaneous return to pre-industrial climate conditions, as surface air temperature and sea level response exhibit a substantial time lag relative to atmospheric CO2

NCIO 1.0: a simple Fortran NetCDF interface

The NetCDF (Network Common Data Form) library has become an indispensable tool for data and model output management in geoscience. However for simple tasks, particularly in Fortran, the complexity of native NetCDF functionality can be cumbersome. The NCIO (NetCDF Input/Output) module has been designed as an interface to the NetCDF library with simplicity and ease of use in mind. While this implies that some NetCDF functionality is masked from the user, the subroutines provided here are adequate for basic serial reading and writing tasks of up to 6-D data arrays along with corresponding data attributes. The code is available online via a GitHub repository (http://www.github.com/alex-robinson/ncio), which includes an example program to illustrate the approach

Simulating the Greenland ice sheet under present-day and palaeo constraints including a new discharge parameterization

© Author(s) 2015. In this paper, we propose a new sub-grid scale parameterization for the ice discharge into the ocean through outlet glaciers and inspect the role of different observational and palaeo constraints for the choice of an optimal set of model parameters. This parameterization was introduced into the polythermal ice-sheet model SICOPOLIS, which is coupled to the regional climate model of intermediate complexity REMBO. Using the coupled model, we performed large ensemble simulations over the last two glacial cycles by varying two major parameters: a melt parameter in the surface melt scheme of REMBO and a discharge scaling parameter in our parameterization of ice discharge. Our empirical constraints are the present-day Greenland ice sheet surface elevation, the surface mass balance partition (ratio between total ice discharge and total precipitation) and the Eemian interglacial elevation drop relative to present day in the vicinity of the NEEM ice core. We show that the ice discharge parameterization enables us to simulate both the correct ice-sheet shape and mass balance partition at the same time without explicitly resolving the Greenland outlet glaciers. For model verification, we compare the simulated total and sectoral ice discharge with other estimates. For the model versions that are consistent with the range of observational and palaeo constraints, our simulated Greenland ice sheet contribution to Eemian sea-level rise relative to present-day amounts to 1.4 m on average (in the range of 0.6 and 2.5 m).Peer Reviewe

Simple models for the simulation of submarine melt for a Greenland glacial system model

Two hundred marine-terminating Greenland outlet glaciers deliver more than half of the annually accumulated ice into the ocean and have played an important role in the Greenland ice sheet mass loss observed since the mid-1990s. Submarine melt may play a crucial role in the mass balance and position of the grounding line of these outlet glaciers. As the ocean warms, it is expected that submarine melt will increase, potentially driving outlet glaciers retreat and contributing to sea level rise. Projections of the future contribution of outlet glaciers to sea level rise are hampered by the necessity to use models with extremely high resolution of the order of a few hundred meters. That requirement in not only demanded when modeling outlet glaciers as a stand alone model but also when coupling them with high-resolution 3-D ocean models. In addition, fjord bathymetry data are mostly missing or inaccurate (errors of several hundreds of meters), which questions the benefit of using computationally expensive 3-D models for future predictions. Here we propose an alternative approach built on the use of a computationally efficient simple model of submarine melt based on turbulent plume theory. We show that such a simple model is in reasonable agreement with several available modeling studies. We performed a suite of experiments to analyze sensitivity of these simple models to model parameters and climate characteristics. We found that the computationally cheap plume model demonstrates qualitatively similar behavior as 3-D general circulation models. To match results of the 3-D models in a quantitative manner, a scaling factor of the order of 1 is needed for the plume models. We applied this approach to model submarine melt for six representative Greenland glaciers and found that the application of a line plume can produce submarine melt compatible with observational data. Our results show that the line plume model is more appropriate than the cone plume model for simulating the average submarine melting of real glaciers in Greenland

Simulating the Greenland ice sheet under present-day and palaeo constraints including a new discharge parameterization

In this paper, we propose a new sub-grid scale parameterization for the ice discharge into the ocean through outlet glaciers and inspect the role of different observational and palaeo constraints for the choice of an optimal set of model parameters. This parameterization was introduced into the polythermal ice-sheet model SICOPOLIS, which is coupled to the regional climate model of intermediate complexity REMBO. Using the coupled model, we performed large ensemble simulations over the last two glacial cycles by varying two major parameters: a melt parameter in the surface melt scheme of REMBO and a discharge scaling parameter in our parameterization of ice discharge. Our empirical constraints are the present-day Greenland ice sheet surface elevation, the surface mass balance partition (ratio between total ice discharge and total precipitation) and the Eemian interglacial elevation drop relative to present day in the vicinity of the NEEM ice core. We show that the ice discharge parameterization enables us to simulate both the correct ice-sheet shape and mass balance partition at the same time without explicitly resolving the Greenland outlet glaciers. For model verification, we compare the simulated total and sectoral ice discharge with other estimates. For the model versions that are consistent with the range of observational and palaeo constraints, our simulated Greenland ice sheet contribution to Eemian sea-level rise relative to present-day amounts to 1.4 m on average (in the range of 0.6 and 2.5 m)